Cementing CO2 into C-S-H: A step toward concrete carbon neutrality.

Addressing the existing gap between currently available mitigation strategies for greenhouse gas emissions associated with ordinary Portland cement production and the 2050 carbon neutrality goal represents a significant challenge. In order to bridge this gap, one potential option is the direct gaseous sequestration and storage of anthropogenic CO2 in concrete through forced carbonate mineralization in both the cementing minerals and their aggregates. To better clarify the potential strategic benefits of these processes, here, we apply an integrated correlative time- and space-resolved Raman microscopy and indentation approach to investigate the underlying mechanisms and chemomechanics of cement carbonation over time scales ranging from the first few hours to several days using bicarbonate-substituted alite as a model system. In these reactions, the carbonation of transient disordered calcium hydroxide particles at the hydration site leads to the formation of a series of calcium carbonate polymorphs including disordered calcium carbonate, ikaite, vaterite, and calcite, which serve as nucleation sites for the formation of a calcium carbonate/calcium-silicate-hydrate (C-S-H) composite, and the subsequent acceleration of the curing process. The results from these studies reveal that in contrast to late-stage cement carbonation processes, these early stage (precure) out-of-equilibrium carbonation reactions do not compromise the material's structural integrity, while allowing significant quantities of CO2 (up to 15 w%) to be incorporated into the cementing matrix. The out-of-equilibrium carbonation of hydrating clinker thus provides an avenue for reducing the environmental footprint of cementitious materials via the uptake and long-term storage of anthropogenic CO2.

Timescale prediction of complex multi-barrier pathways using flux sampling molecular dynamics and 1D kinetic integration: Application to cellulose dehydration.

Reactive molecular dynamics (MD) simulations, especially those employing acceleration techniques, can provide useful insights on the mechanism underlying the transformation of buried organic matter, yet, so far, it remains extremely difficult to predict the time scales associated with these processes at moderate temperatures (i.e., when such time scales are considerably larger than those accessible to MD). We propose here an accelerated method based on flux sampling and kinetic integration along a 1D order parameter that can considerably extend the accessible time scales. We demonstrate the utility of this technique in an application to the dehydration of crystalline cellulose at temperatures ranging from 1900 K to 1500 K. The full decomposition is obtained at all temperatures apart from T = 1500 K, showing the same distribution of the main volatiles (H2O, CO, and CO2) as recently obtained using replica exchange molecular dynamics. The kinetics of the process is well fitted with an Arrhenius law with Ea = 93 kcal/mol and k0 = 9 × 1019 s-1, which are somehow larger than experimental reports. Unexpectedly, the process seems to considerably slow down at lower temperatures, severely departing from the Arrhenius regime, probably because of an inadequate choice of the order parameter. Nevertheless, we show that the proposed method allows considerable time sampling at low temperatures compared to conventional MD.

Poroelasticity of Methane-Loaded Mature and Immature Kerogen from Molecular Simulations.

While hydrocarbon expulsion from kerogen is certainly the key step in shale oil/gas recovery, the poromechanical couplings governing this desorption process, taking place under a significant pressure gradient, are still poorly understood. Especially, most molecular simulation investigations of hydrocarbon adsorption and transport in kerogen have so far been performed under the rigid matrix approximation, implying that the pore space is independent of pressure, temperature, and fluid loading, or in other words, neglecting poromechanics. Here, using two hydrogenated porous carbon models as proxies for immature and overmature kerogen, that is, highly aliphatic hydrogen-rich vs highly aromatic hydrogen-poor models, we perform an extensive molecular-dynamics-based investigation of the evolution of the poroelastic properties of those matrices with respect to temperature, external pressure, and methane loading as a prototype alkane molecule. The rigid matrix approximation is shown to hold reasonably well for overmature kerogen even though accounting for flexibility has allowed us to observe the well-known small volume contraction at low fluid loading and temperature. Our results demonstrate that immature kerogen is highly deformable. Within the ranges of conditions considered in this work, its density can double and its accessible porosity (to a methane molecule) can increase from 0 to ∼30%. We also show that these deformations are significantly nonaffine (i.e., nonhomogeneous), especially upon fluid adsorption or desorption.

From cellulose to kerogen: molecular simulation of a geological process.

The process by which organic matter decomposes deep underground to form petroleum and its underlying kerogen matrix has so far remained a no man's land to theoreticians, largely because of the geological (Myears) timescale associated with the process. Using reactive molecular dynamics and an accelerated simulation framework, the replica exchange molecular dynamics method, we simulate the full transformation of cellulose into kerogen and its associated fluid phase under prevailing geological conditions. We observe in sequence the fragmentation of the cellulose crystal and production of water, the development of an unsaturated aliphatic macromolecular phase and its aromatization. The composition of the solid residue along the maturation pathway strictly follows what is observed for natural type III kerogen and for artificially matured samples under confined conditions. After expulsion of the fluid phase, the obtained microporous kerogen possesses the structure, texture, density, porosity and stiffness observed for mature type III kerogen and a microporous carbon obtained by saccharose pyrolysis at low temperature. As expected for this variety of precursor, the main resulting hydrocarbon is methane. The present work thus demonstrates that molecular simulations can now be used to assess, almost quantitatively, such complex chemical processes as petrogenesis in fossil reservoirs and, more generally, the possible conversion of any natural product into bio-sourced materials and/or fuel.

Cement As a Waste Form for Nuclear Fission Products: The Case of (90)Sr and Its Daughters.

One of the main challenges faced by the nuclear industry is the long-term confinement of nuclear waste. Because it is inexpensive and easy to manufacture, cement is the material of choice to store large volumes of radioactive materials, in particular the low-level medium-lived fission products. It is therefore of utmost importance to assess the chemical and structural stability of cement containing radioactive species. Here, we use ab initio calculations based on density functional theory (DFT) to study the effects of (90)Sr insertion and decay in C-S-H (calcium-silicate-hydrate) in order to test the ability of cement to trap and hold this radioactive fission product and to investigate the consequences of its ß-decay on the cement paste structure. We show that (90)Sr is stable when it substitutes the Ca(2+) ions in C-S-H, and so is its daughter nucleus (90)Y after ß-decay. Interestingly, (90)Zr, daughter of (90)Y and final product in the decay sequence, is found to be unstable compared to the bulk phase of the element at zero K but stable when compared to the solvated ion in water. Therefore, cement appears as a suitable waste form for (90)Sr storage.